Asymmetric Berry-Phase Interference Patterns in a Single-Molecule Magnet 
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A Mn4 single-molecule magnet displays asymmetric Berry-phase interference patterns in the 
transverse- field {Ht) dependence of the magnetization tunneling probability when a longitudinal 
field {Hl) is present, contrary to symmetric patterns observed for Hl — 0. Reversal of Hl results 
in a reflection of the transverse-field asymmetry about Ht = 0, as expected on the basis of the 
time-reversal invariance of the spin-orbit Hamiltonian which is responsible for the tunneling oscilla- 
tions. A fascinating motion of Berry-phase minima within the transverse field magnitude-direction 
phase space results from a competition between non-collinear magneto- anisot ropy tensors at the two 
distinct Mn sites. 

PACS numbers: 75.45.+j, 75.50.Xx 



Almost two decades of research have established single- 
molecule magnets (SMMs) as prototype systems for un- 
derstanding fundamental quantum phenomena associ- 
ated with nanoscale magnetism [H, Q , as well as demon- 
strating their potential for future applications [3]. The 
most important characteristics of SMMs can be modeled 
reliably using a giant-spin approximation (GSA) whereby 
the molecule is treated as a rigid magnetic unit with total 
spin, 5, weakly interacting with its enviroment. Indeed, 
this model accounts for the essential features of the quan- 
tum tunneling of magnetization (QTM) observed in these 
molecules, as well as its quenching due to Berry phase 
interference (BPI) resulting from different tunneling tra- 
jectories [j-l6|. However, the GSA ignores the internal 
couplings within a SMM, thereby completely failing to 
account for QTM transitions that involve fluctuations of 
the total spin of the molecule 0-[2| , or otherwise obscur- 
ing intrinsic relationships that exist between QTM se- 
lection rules and the underlying molecular structure [lo| . 
In a recent interesting example, it was demonstrated that 
a tilting of the zero-field splitting (zfs) tensors in a tri- 
angular Mug^^ SMM (lowering the symmetry of the spin 
Hamiltonian from Cq to C3) results in new QTM selec- 
tion rules and strongly affects the transverse field depen- 
dence of the remaining forbidden QTM resonances |1Q| . 
These observations likely explain the absence of QTM se- 
lection rules in most SMMs studied to date, since internal 
dipolar fields and/or weak sample disorder are often suf- 
ficient to cause observable relaxation at resonances oth- 
erwise forbidden by symmetry. A detailed understanding 
of these and related phenomena has mostly been facili- 
tated by studies of low nuclearity SMMs [ill, |T2| , where 
exact diagonalization of the multi-spin Hamiltonian en- 
ables consideration of the internal degrees of freedom of 
the molecule. 

In this letter, we focus on the QTM relaxation asso- 
ciated with a centro-symmetric mixed- valent Mn2^Mn2^^ 




FIG. 1: (Color on-line) Hysteresis loops recorded as a function 
of Hl at different temperatures. The inset shows the relevant 
magnetic centers (i.e. Mn^^, Mn^^^, N and O) at the core of the 
Mn4-Bet unit. The principal magnetic axes of the molecule 
are indicated, with the easy {z-) axis determined by the JT 
elongation associated with the Mn ions (Mn"^-N line, see 
main text); the ^-axis lies in the plane defined by the two 
Mn^^^ ions and their JT axes. 

complex which shows an asymmetric BPI pattern with 
respect to the polarity of the transverse component of the 
applied field {Ht ^ magnetic easy axis). We show that 
this behavior results from a competition between non- 
collinear magneto-anisotropy tensors at the two crystal- 
lographically distinct Mn ions, which is also responsible 
for an unusual motion of the Berry-phase minima within 
the transverse field magnitude-direction phase space. We 
show how the asymmetry can be inverted upon reversal 
of the longitudinal field {Hl \\ easy axis), i.e., the BPI 
pattern is invariant with respect to a full inversion of the 
applied field, consistent with the time-reversal symmetry 
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of the underlying zero-field Hamiltonian. 

The [Mn4 (Bet)4 (mdea)2 (Hmdea)2] (BPh4)4 complex 
(henceforth Mn4-Bet) crystallizes in the triclinic PI 
space group with half the molecule in the asymmetric 
unit; the other half is generated via inversion, result- 
ing in the four Mn ions lying in a plane (the molecular 
plane), with the Mn^^^ Jahn- Teller (JT) axes oriented 
along the Mn^^^-N bonds, which lie 122.61 degrees out of 
this plane, i.e., roughly perpendicular to the molecular 



plane [13|, [l^. There are no solvent groups in the lat- 



tice and the four BPh4 anions enhance isolation, result- 
ing in extremely clean X-ray diffraction and EPR data 
14, lH. A sketch of the Mn4 core, where the magnetic 
axes are indicated, is inset to Fig. 1. Magnetic and EPR 
measurements performed at relatively high temperatures 
(T > 2 K) suggest a spin S = 9 ground state, and that 
Mn4-Bet is a SMM with a barrier of - 20 K 0. 

A high-sensititivity micro-Hall effect magnetometer, a 
He^/He^ dilution fridge and a 3D vector superconducting 
magnet were employed to record magnetization hysteresis 
curves as a function of a magnetic field applied parallel 
to the easy axis of the molecules [l6|, at temperatures 
down to 35 mK. The results are shown in Fig. 1, where 
extremely sharp QTM resonances (/c = 0, 1 & 2), spaced 
by AHl ^ 0.21 T, confirm the high quality of the crystal. 
Within the GSA, this spacing corresponds to an axial zfs 
parameter, D = —0.28 K (^ = 2). The observed blocking 
and crossover temperatures are ~1.2 K and ~0.2 K, re- 
spectively. A transverse field was subsequently employed 
in order to study the symmetry of the QTM in resonances 
k = and k = 1. Fig. 2a shows the modulation of the 
QTM probability, = {Mf - Mi)/{Msat - Mi) [13], 
for resonance /c = 0, as a function of Ht applied along 
the magnetic hard axis {(f) = 0°). This angle, which lies 
^ 30° away from one of the crystal faces, was deduced 
from the two-fold modulation of Pq as a function of the 
orientation, ^, of a 0.2 T transverse field within the hard 
plane (see inset to Fig. 2a) [Tij . 

The Pq oscillations in Fig. 2a correspond to BPI, with 
minima at regularly spaced field values {AHt = 0.3T). 
A maximum in Pq is found at Ht = 0, as expected 
for an integer spin value. Within the GSA, AHt = 
2kB{2E[E^D]f^'^/giiB [iQ^, yielding a 2^^-order rhom- 
bic zfs parameter, E = ±60 mK. Note that the regularly 
spaced k = BPI minima are invariant under inversion 
of Ht^ i.e., they are symmetric with respect to Ht = 0. 
Interestingly, this is not the case in resonance /c = 1, 
for which the behavior of the QTM probability is very 
different. This can be seen in Fig. 2b, which illustrates 
the dependence of Pi on Ht^ for (j) = 13.5° (the angle 
for which the first BPI minimum at Ht = 0.30 T is the 
sharpest). In fact, for resonance /c = 1, different BPI 
minima appear at different field orientations, (/), of the 
transverse field within the xy (hard) plane of the molecule 
(see Fig. 3) fl^, i.e., the first minimum {Ht = 0.3 T) ap- 
pears at = 13.5°, while the second {Ht = 0.6 T) occurs 
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FIG. 2: (Color on-line) Modulation of the QTM probabilities 
for resonances /c = (a) and /c = 1 (b) as a function of Ht 
applied at different angles, 0, within the hard plane of the 
Mn4-Bet SMM 0. The inset to (a) illustrates the two-fold 
angular modulation of Pq for Ht = 0.2 T, providing clear 
evidence for a significant 2^^-order rhombic anisotropy. The 
asymmetry of the BPI pattern of oscillations in resonance 
/c = 1 is inverted upon reversal of Hl- The inset to (b) 
illustrates the classical anisotropy barrier generated by the 
non-collinear zfs tensors (see main text for explanation) and 
the different perspectives resulting from permutations of ±Ht 
and z^Hl- 



at (j) = 6°, contrary to what is found for the k = res- 
onance (all Pq minima are seen most clearly at ^ = 0°) . 
Such behavior has been predicted theoretically 2qI. [21|. 
though never observed experimentally. 

Before considering this aspect in detail, we first discuss 
the asymmetric nature of the BPI oscillation pattern in 
resonance k = 1. As seen clearly in Fig. 2b, reversal of 
the longitudinal field, Hl^ results in a reflection of the 
Pi BPI pattern about Ht = 0. In other words, the BPI 
minima are in fact invariant under a full magnetic field 
inversion, as required on the basis of the time-reversal 
invariance of the spin-orbit Hamiltonian responsible for 
this physics. As noted above, the symmetries of BPI pat- 
terns must respect the symmetry of the zero-field spin 
Hamiltonian. If one considers only 2"^-order zfs within 
the GSA, then the resulting Hamiltonian necessarily be- 
longs to the orthorhombic point group and possesses the 
following symmetry elements: (1) three mutually orthog- 
onal two- fold rotation axes (x, y h z)] (2) three mutually 
orthogonal mirror planes {xy^ xz & yz); and (3) an in- 
version center. (2) guarantees invariance with respect 
to reversal of Ht, i.e., it enforces symmetric BPI pat- 
terns, irrespective of whether a longitudinal field is ap- 
plied {k > 0) or not {k = 0). As we show below, one must 
break the xy mirror symmetry in order to obtain asym- 
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metric BPI patterns with respect to inversion of Ht- In 
this case, reversal of Hl results in different patterns; the 
time reversal symmetry then guarantees that these two 
patterns are mirror images. Nevertheless, no matter how 
many spatial symmetries are broken, the time-reversal 
invariance of the spin-orbit interaction guarantees that 
the BPI minima should be invariant under a full reversal 
of the applied field, i.e., simultaneous reversal Hl and 
Ht-, as we have confirmed experimentally. 

It is possible to reproduce the essential features of 
the experiments by introducing 4*^-order terms into 
the GSA; the xy mirror symmetry can then be broken 
by rotating the coordinate frames of the 2"^ and 4*^- 
order tensors. Interestingly, this approach also repro- 
duces the complex motion of the Pi minima within the 
Ht — 4> phase-space shown in Fig 3. However, a com- 
plete GSA analysis for Mn4-Bet requires many param- 
eters and provides little insight, while the same physics 
can be naturally understood within a multispin descrip- 
tion which considers the internal structure of the Mn4- 
Bet molecule. Note that the emergence of significant 
higher-order anisotropy terms within the GSA is a mani- 
festation of mixing of the ground spin state with excited 
states, which can only be captured within a multi-spin 



model [UJ, ll2|. In this context, the xy mirror sym- 



metry may be trivially broken by rotating (tilting) the 
zfs tensors at the two inequivalent magnetic sites in the 
molecule so that their local 2;-axes no longer coincide (this 
is similar to the case of the Mug^^ SMM discussed in the 
introduction [TqI,!!!). 

To explain the experimental findings we have diago- 
nalized the multi-spin Hamiltonian, where the four Mn 
ions are coupled according to the sketch in Fig. 4a. 

3 

H = ^{si • R[ ' di' Ri'Si- giiBSi • B) 



s,'A-Sj). (1) 



The first term represents the local anisotropy of the i^^ 
ion, di being the 2"^-order zfs tensor {di^xx = e^, di^yy = 
—Ci & di^zz = di^ with di & representing the axial and 
rhombic anisotropics, respectively); the Ri are Euler ma- 
trices (defined by Euler rotation angles ai^ j3i h ^i) spec- 
ifying the orientations of these tensors. The second term 
is the Zeeman coupling to the applied field, where we as- 
sume an isotropic Lande factor, g = 2.00. The third and 
fourth terms represent the exchange and dipolar inter- 
actions, respectively. These terms are also time-reversal 
invariant, and do not change any of the preceeding argu- 
ments. As depicted in Fig. 4a, three independent near- 
neighbor exchange coupling constants, Ja, Jh and Jc, are 
considered: Ji,2 = ^3,4 = Ja-, -^2,3 = Ji,a = Jh-, -^2,4 = Jc 
and Ji^3 = 0. The dipolar matrix, j has been chosen 
to exactly reproduce all of the dipolar couplings within 
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FIG. 3: (Color on-line) Contour plots of the QTM probabil- 
ities for resonances k — ^ and /c = 1 as a function of Ht 
and (j). All of the k — minima lie approximately along the 
= 0° axis, whereas the k — 1 minima appear at different 
orientations for different Ht values. 



the molecule, with no fitting parameters. 

Figure 4b plots the locations of BPI minima obtained 
via diagonalization of Eqn. [1] (solid red symbols). This 
simulation, which takes into account the small misaligne- 
ment of the experimental field rotation plane [l8|, em- 
ployed the following parameters: ^2 = <^4 = —4.99 K and 
62 = 64 = 0.82 K, with the easy and hard anisotropy 
axes along z (a2 = 0) and x {P2 = 0), respectively; 
di = ds = —0.67 K & ei = 63 = 0, with the axes rotated 
with respect to the central spin by identical Euler angles 
<^i,3 = 45°, Pi^s = 0° (as required by inversion symme- 
try); 7 being zero for all ions; finally, isotropic ferromag- 
netic exchange constants Ja = —3.84 K, J5 = —1.20 K 
and Jc = —3.36 K are used. It should be stressed that 
these parameters are additionally constrained by the lo- 
cations of hysteresis loop steps (Fig. 1) and extensive 
angle-dependent EPR measurements [l5|. Moreover, the 
obtained anisotropy values for the Mn^^^ ions are very 
similar to related Mn^^^ complexes [10], while the di^s 
value lies within the bounds reported for other Mn^^ sys- 
tems [23I. The quantitative agreement with experiment 
is also excellent. The motion of the Pi minima can be 
understood as a result of the competition between dif- 
ferent anisotropic interactions within the molecule, with- 
out a need to invoke unphysical 4*^ and higher order 
anisotropics. Importantly, the angular positions (0) of 
the k = 1 minima move with Ht^ while the k = min- 
ima remain stationary, as found experimentally (Fig. 3). 

Finally, the multi-spin model perfectly reproduces the 
Ht asymmetry of the k = 1 BPI pattern. As seen in 
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FIG. 4: (Color on-line) (a) Sketch of the Mn4 core indicat- 
ing the interaction parameters used to explain the results, 
(b) Measured (open black symbols) and calculated (solid red 
symbols) ^/i^T-dependence of the BPI minima for resonances 
A; = (circles) and 1 (squares) obtained from Eqn.[T] (c) Cal- 
culated tunnel splittings for resonance k = 1, for = 9.5°, 
as a function of Ht for Hl > (continuous balck line) and 
Hl < (discontinuous blue line). 

Fig. 4c, the asymmetry is reversed upon inversion of 
Hl^ as required by the time-reversal invariance of the 
anisotropic interactions in Eqn. (TJ and observed experi- 
mentally (Fig. 2b). The crucial ingredient is the tilting 
of the zfs tensors of the external spins, si & 53, relative to 
the central spins 52 & 54, so that the xy mirror symme- 
try is broken. This is illustrated in the inset to Fig. 2b, 
where one observes that the classical energy landscape 
is invariant under full field inversion (blue vs. black ar- 
rows), while this is not the case when only Ht is reversed 
(red vs. black arrows). The Euler angle ai^s = 45° 
results in a significant projection of the relatively weak 
anisotropy associated with the Mn^^ ions into the hard 
{xy) plane. This, together with the finite 62,4 param- 
eters and the dipolar interactions, results in competing 
transverse interactions and to the complexity of the BPI 
patterns observed in Fig. 3. We note that the dipolar in- 
teraction has a very significant effect on the energy levels 
of the molecule: the zfs within the S = 9 multiplet varies 
by as much as 10% when dipolar interactions are omit- 
ted, and the location of the k = 1 QTM step is shifted 
by - 0.02T. 

We conclude by noting that asymmetric BPI patterns 
have been seen in other centro-symmetric SMMs for 
which a clear explanation has been lacking 0, [sl, H^. 
The present results may help shed light on the effect that 
symmetric anisotropic interactions can have in magnetic 
systems with inversion symmetry, where a net antisym- 
metric interaction is strickly forbidden. The present work 
clearly demonstrates how studies of simple low nuclear- 
ity systems can address fundamental symmetry consid- 



erations related to QTM in molecular nanomagnetism. 
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